NMR relaxation rates and Knight shifts in MgB 2 
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We calculate ab initio the NMR relaxation rates and the Knight shifts in MgB2. We show 
that the dominant relaxation mechanism at the n B nucleus is the interaction with the electronic 
orbital moment, and we give a simple explanation of that using a simple sp tight binding model. 
When Stoner enhancement (also calculated ab-initio) is accounted for, we obtain good agreement 
with reported experimental values. For the 25 Mg nucleus, we predict that the dominant relaxation 
mechanism is the Fermi-contact interaction, which also dominates the Mg Knight shift. 
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Recent discovery^! of superconductivity in MgB2 cre- 
ated substantial interest. It was suggested that the un- 
derlying mechanism is electron-phonon interaction in the 
boron sublatticeo, which was subsequently confirmed by. 
observation of a sizeable boronoo, but not magnesiumcl 
isotope effect. State of the^arf local density approxi- 
mation (LDA) calculationsuLJ produced electron-phonon 
coupling constants A ranging from 0.75 to 0.87. Lacking 
single crystals, experimental determination of Arselies on 
the specific heat renormalization measurementstl. Using 
the LDA density of states (DOS), these experiments give 
A ~ 0.6 — 0.8; however, if there is any many-body renor- 
malization of the LDA DOS, these experiments should 
be reanalyzed. 

Nuclear magnetic resonance (NMR) is a common probe 
of the DOS. The measured quantities, the spin-lattice 
relaxation rate, 1/Ti, and the Knight shift, K, are re- 
lated to the spin susceptibility, and thus are not sub- 
ject to a phonon renormalization. Measurements of the 
relaxation rates and the Knight shift of 11 B already 
existH E3. FroHt-the electronic structure of MgB2 one 
can conjectureou3 that the main source of relaxation 
should be the hypcrfine coupling between the nuclear 
spin and conduction B p electrons. However, a full mi- 
croscopic understanding of the NMR data (relaxation 
rates and_Knight shifts) is still missing. While different 
sourcesO'ErEj reasonably agree among themselves about 
the relaxation rates, reporting 1/T%T between 5.6 xlO -3 
and 6.5xl0~ 3 1/(K sec), there is considerabLe,contro- 
versy about the Knight shifts. Some authorsEd report 
a small average shift K = (K z + 2K xy )/3 = 0.0175%, 
and give an upper bound on its anisoteppy, K ax = 
(K z - K xy )/3 < 0.0030%. Other authors^ report even 
smaller (K = 0.006%) shift and they attribute the shift 
to the Fermi-contact interaction. Note that the Korringa 
relation, r = K 2 (T 1 T)(j n /2n B ) 2 (47rk B h) « 1, where j n 
is the nuclear gyromagnetic ratio, is not satisfied here, 
as the measurements give r ~ 0.2. Finally, a tiny neg- 
ative shift (K — —0.0005%) was measured by Tou et 
al.Jlj and attributed to core polarization. These discrep- 



ancies might arise from the difficulties in measuring the 
n B shift, due to its smallness, andHpossibly, from the se- 
lection of the reference material £rfc3 Therefore, in order 
to clarify the microscopic origin of the NMR relaxation 
process and of the Knight shift, ab-initio calculations are 
highly desirable. 

In the present work we report LDA calculation of the 
relaxation rates and of the Knight shifts. We will show 
that for n B the relaxation is due to the p states, and the 
orbital relaxation rate is about 3 times larger than the 
dipole rate and 10 time larger than the Fermi- contact 
rate. After an appropriate Stoner renormalization is in- 
cluded, the agreement with the experiment is very good. 
On the other hand, the main source of Knight shift is the 
hyperfine coupling with s electrons. Also, the (yet un- 
measured) relaxation on Mg is mainly due to the Fermi- 
contact interaction with the s states. 

The hyperfine interaction — ^7 ?l I • H is the coupling 
between the nuclear magnetic moment H'y n l and the hy- 
perfine field H produced at the site of the nucleus by 
the conduction electrons. In order to discuss separately 
the different relaxation mechanisms, we neglect the small 
spin orbit coupling and split the hyperfine interaction 
into three terms, -h r y n L[H. +H d +H F }. The first term is 
the coupling with the electronic orbital moment; the sec- 
ond and the third terms are, respectively, the dipole and 
the Fermi-contact interaction with the electronic spin. 
Thus the total hyperfine field is given by 
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where r, s and 1 are the electronic position, spin, and an- 
gular momentum operator. In the case of 11 B, / = 3/2 
and 7„ = 0.89 7jv, while in the case of 25 Mg 1 — 5/2 and 
7„ = — 0.177AT, with 7tv = e/m p c. 

According to Fermi's .golden rule, the relaxation rate, 
1/Ti, rnay be written aa!3 
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where /(e) is the Fermi-Dirac distribution, s is the spin 
index, and |m) are the eigenstates of I z . Here k stands 
for both the wave vector and band index. Expansion 
of the Fermi function and integration over the nuclear 
spin yields^ for a polycrystalline sample, the following 
expression^ 
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where N is understood as a diagonal matrix in the 
spin and k-space, 8 ss '0~)sk'&{tv)- The relevant prefac- 
tor is C — (47r/cs/fi,) (7„/7 e ) 2 , the same that appears 
in the Korringa relation. In the present case C ~ 
1.4 x 10 4 /(K sec) for 25 Mg and C ~ 3.9 x 10 5 /(K sec) 
for 11 B. The interaction cross terms in Eq. ([!]), i.e. the 
terms proportional to Tr[H°7VH d 7V], Tr[H°iVH F 7V] and 
Tr[H F iVH d iV] all vanish, the first two exactly, because 
Tr[s]=0 and the third vanishes exactly for polycrystals 
because Tr[s 2 — 3(s • r) 2 ] = 0, and approximately for sin- 
gle crystals, when the d-electron DOS is small(c/. Ref. 
Eq , Eq. 24). Thus, without the core polarization, which 
will be discussed later, the relaxation rate has three con- 
tributions: the orbital, the dipolc and the contact-field 
term. Note that in the terminology of Ref.Ej, all cross 
terms, diagonal in interaction but off-diagonal in angular 
momentum, are included in the calculation. More details 
on this derivation can be found in Ref. [l^. 

In order to evaluate the relaxation rate, we adopt the 
tight buiding LMTO-ASA method (LMT047 Stuttgart 
code)o This method has been already used with suc- 
cess to calculate 1/Ti, e.g. in AaCeo-tJ Thus we express 



the Bloch function as |iks) 
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4>Ri{ti>RU t)Yl(tr), where 4>Ri is the radial solution of the 
Schrodinger equation at the energy c^rl, 4>rv is its en- 
ergy derivative, Yl is a spherical harmonic with L = Im. 
For simplicity, in the following we will write only the con- 
tributions from 4>ri, although in the calculation we have, 
of course, included all terms. Thus the three contribu- 
tions to 1/Ti can be expressed as 
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and of the radial integrals involving 4>Ri(£^Ri,r) 



(r 3 )m = J 4>Ri{e vR i,r)r 3 (j)Rv{e V Ri' ,r)r 2 dr. (3) 

The Fermi-contact, the orbital, and the dipole contribu- 
tions may then be written respectively as 
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Here D» L , = (L%\L), h = h, l±i 
Clu = ^fY 2fl (r)Y L (r)*Y L ,(r)£r 

In the same way, the Knight shift can be written as 
K a =2fiB Tr(1 IHqA^I f) where a is the direction along 
which the external magnetic field is applied. As the re- 
laxation rate, the relative shift may also be expressed as 
a function of the DOS matrix and the radial integrals, 
expanding the Bloch function in the LMTO basis set. 

The DOS matrix was calculated by the linear tetrahe- 
dron method. We found that the results were already 
very well converged with a mesh of 370 irreducible k 
points. In order to minimize the linearization error ac- 
curate wavefunctions at the Fermi level, the final runs 
were performed with t V Ri = ep. The convergence of the 
sums over the angular momentum was also very good. 
We find that we can truncate after 1 = 2. The reason is 
that the radial integrals ((ao/r) 3 )/;/ (ao is the Bohr ra- 
dius) decrease quickly when I and I' increases. For Mg 
we find, e.g. ((a /r) 3 )n = 4.8, ((a /r) 3 ) 22 = 0.16, and 
((aoA) 3 >33 = 0.09, and in for B, ((a /r) 3 )n = 1.1 and 
((aoA) 3 > 22 = 0.2. 

What is the dominant mechanism that gives rise to 
the magnetic relaxation at B and Mg nuclei? In most 
metals it is the Fermi contact one, defined by the DOS 
of the s electrons at the Fermi level. However, in the 
case of MgB 2 the states near the Fermi level are mainly 
B p. We find that the ratio N ss (Mg)/N tot (Mg) - 1/4, 
and N ss (B)/N tot (B) ~ 1/50. Therefore, at least in the 
case of B, the ratio is very small, and the Fermi contact 
term could become comparable or even smaller than the 
dipole or the orbital term. We have calculated all three 
contributions for both elements and show the results in 
Tables 1 and 2. 

We also calculate ab initio the core polarization. For 
this purpose we applied in the calculations an exter- 
nal magnetic field B, and then calculated m„(0), the 
spin density of the nth core shell at the nucleus. Then 
the core polarization Knight shift can be obtained as 
K cp = Ms(8tt/3) X)„( m "(0)/-B), and the corresponding 
contributions to the relaxations rate can be computed 
from the Korringa relations for the core statesllS. For 
25 Mg, the contribution of the core polarization is negli- 
gible and the Fermi-contact interaction dominates. For 
11 B, we find that the contribution of the Is shell is of the 
same order, but of opposite sign, as the the 2s shell con- 
tribution. Their total effect is thus small. Also, since the 
Fermi-contact contribution for B is much smaller than in 
Mg, the relative effect of the dipole term is larger, lead- 
ing to a noticeable anisotropy of the Knight shift (about 
30%), while the Mg Knight shift is essentially isotropic. 
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TABLE I. Knight shift, Ka in %. Both unrenormalized and Stoner-enhanced values are included, as discussed in the text. 
The label a =ab,c indicates the direction of the external magnetic field. 





dipole (ab) 


dipole (c) 


orbital 


Fermi- contact 


core 


Total (ab/c) 


Total (renormalized) 


Expt. a 


Expt. b Expt. c 


Mg 


0.0005 


-O.OOfO 





0.0260 


0.0003 


0.027f/0.0256 


0.0361/0.0341 






B 


-0.0004 


0.0008 





0.0027 


-0.0007 


0.0016/0.0028 


0.0024/0.0042 


0.0175 


0.006 -0.0005 



''Ref. b Ref. § c Ref. ^ 



TABLE II. Relaxation rate 1/TiT in 10~ 3 /( K sec). Both unrenormalized and Stoner-enhanced values are included, as 
discussed in the text. 





orbital 


dipole 


Fermi-contact 


core 


Total 


Total (renormalized) 


Expt. a 


Expt. b 


Expt. c 


Mg 


0.02 


0.01 


1.0 


0.0001 


1.0 


1.6 








B 


2.6 


0.8 


0.28 


0.02 


3.7 


5.9 


5.6 


6.5 


6.1 



"Ref. ^ b Ref. [H]; c Ref. | 



In order to understand the numerical results, we first 
calculate analytically the shifts and the relaxation rate 
for a model Hamiltonian which includes only B s and B p 
orbitals. We start with the contribution of s electrons, i.e. 
the contact term. The contact shift may be written as 
K ~ p 2 B (A/3)\4> s (0)\ 2 N ss , where N ss is the s-projected 
DOS per atom per spin. We find |</> s (0)| 2 /(47r) — 
2.16ap 3 , and N ss ~ 0.002 states/eV per B atom. There- 
fore K ~ 0.0026% and 1/TiT - 0.25 x 10" 3 1/(K sec). 
Both numbers are very close to those obtained from the 
full calculations (Tables 1,2). 

We now consider the contribution of B p electrons. 
The states at the Fermi level are ~ 70% 
N„_ „_ ~ 0.035 states/(spin eV atom), and 



B 

N, 



p-like. 



0.045 states/(spin eV atom). Thus N PxtPt 
N p /3, where N p is the total p-projected DOS per spin 
per atom. Therefore we find the following approximate 
expression of the orbital contribution to the relaxation 
rate 
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where Tr m [1 • 1] = (21+1)1(1+1) = 6. We find ((a /r) 3 ) ~ 
1.14 and N p /3 ~ 0.038 states/eV per B atom, and there- 
fore 1/TiT - 3 x 10~ 3 /(K sec). The orbital part of the 
Knight shift is zero in this model because nondiagonal 
elements of the DOS matrix vanish. 

In the same way the dipole term can be written as 
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where 



A C lm,lm') 2 = 6 / 5 - Thus we find that the 
B p electron contribution to the dipole relaxation rate 

is 1/TxT - 0.8 x 10" 3 1/(K sec). For the Knight shift 

we find K* 2K T 



, xy - 2 M | C? %(r- 3 ) n (AW3) ~ 
0.0015%. Again, all these numbers are rather close to 
the all-electron results shown in the Tables. The ratio 

(TiW(Ti) or6 ~ (2/3)Tr[l-l]/E /im ,„ 



3.3. The reason for which the orbital term dominates 
over the dipole term is that all three p are present at the 
Fermi level, as opposed, for instance, to the fullerenesEj, 
where only one orbital is present at the Fermi level and 
thus the orbital moment is quenched. 

In terms of the linear response theory, both the 
Knight shift and the relaxation- rate are defined by 
the electronic spin susceptibilitytll, %(q, cj), specifically, 
K cx Rex(O,0), and 1/Ti oc lim^ E q Imx(q, 
Electron-hole excitations renormalized the spin suscep- 
tibility, and in the simplest possible approximation one 
writes 

X(q,w) «xo(q,^)/[l - Jxo(q,w)], 

where xo is the bare (noninteracting) susceptibility, / is 
the so-called Stoner factor, characterizing intraatomic ex- 
change, and the calculations described above correspond 
to total neglect of the Stoner renormalization. One can 
estimate / from LSDA calculation with fixed total spin 
moment by fitting the total energy to the Stoner expres- 
sion, E tot (M) = M 2 /4:N - M 2 I/4: 7 where M is the spin 
moment and N is the total DOS per spin. In this way, 
we found IN = Ix(®> 0) ~ 0.25. Thus we can estimate 
renormalized Knight shift as K sa Kq/(1— IN) — 1.33iio- 
The renormalized values are also shown in the Tables. 

Renormalization of 1 /Tt is somewhat more difficult to 
take into account. It is easy to showc3 that in the Stoner 
approximation 



Imx(q,w) w!mXo(q,w)/[l - 7Rexo(q,w)]' 
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however, averaging this expression over q's requires 
knowledge of the q-dependence of xo- Generally speak- 
ing, renormalization factor lies between 1/(1 — IN) and 
1/(1 — IN) 2 . Using the Lindhard susceptibility ,andr-a 
sphere for the Fermi surface, Shastry and Abrahamsc2l 
found that in the 3D case 



Imxo(q,^ 
i [1 - JRexo(q, 




(Imxo(q,^)) 
(1 - IN)(1 - 2JN/3)' 
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which is a good approximation for IN < 0.7. By 
integrating numerically Eqj7| with the Lindhard func- 
tion, we found a better approximation, good for essen- 
tially all IN, and preserving the correct small IN limit, 
namely (Imxo(q, w )) /(I ~ I-/V) 5 / 3 .Thus we used the fac- 
tor 1.33 5/3 wl.6 for l/T x . For the 2D free electron gas, 
there is no q-dependence in for q < 2fcp, and 

thus the renormalization factor is 1/(1 — IN) 2 . 

The value of 1/T X T = 5.9 x 1(T 3 1/(K sec), is in a 
good agreement with the reported experimental num- 
ber. This means that the DOS, calculated within LDA, 
is a good approximation (maybe a slight underestimate) 
of the bare DOS, and thus the values for the electron- 
phonon coupling constant A, obtained from the specific 
heat measurements, are reliable. 

To the best of our knowledge, there are at present no 
experimental data for Mg. We predict that the mag- 
netic shift is isotropic and that the principal relaxation 
mechanism is the Fermi-contact interaction, despite of 
the fact that N ss /J2i>o^ u ~ V^- The reason is that 
the quantities that one has to compare are not the partial 
DOS N ss and Nu but rather the dimensionless couplings 



(2^/3)1^(0)1^ and M |£ i>0 < 



iNu, and thus the 



relevant ratio is R = (2/3)\MWN ss /(J2i >0 (^ 3 )uNii). 
In the case of Mg we find |0 s (O)| 2 /4tt = 4.54a~ 3 , and 
(f~ 3 )n = 4.8<2q~ 3 . Hence we find R ~ 5. Instead, in the 
case of B, |0 s (O)| 2 /4tt = 2.16a 3 and (r~ 3 )ii = l.la^ 3 , 
and thus R ~ 0.35. The coupling with non s elec- 
tron competes with or dominates over the coupling with 
selectrons when R < 1 

Finally, we would like to mention that the presented 
values for 1/T\T include contributions from both quasi- 
2D p a and 3D p^ bands. If, as suggested!], two different 
gaps open below T c in these bands, the temperature de- 
pendence of \jT\T at low temperature should be com- 
puted taking the different character of these bands in 
the normal states. It is not obvious apriori that the cor- 
responding weights will be just the densities of states. 
Calculations similar to those described above, but band- 
decomposed are needed. 

To summarize, we report first-principles calculations of 
the NMR relaxation rates and the Knight shifts on both 
sites in MgB2. The results are in a good agreement with 
the experiment, provided that the dipole and the orbital 
hyperfine interactions are taken into account, as well as 
the Stoner renormalization of susceptibility. NMR relax- 
ation at n B nucleus is dominated by the orbital interac- 
tion, and that at the 25 Mg nucleus by the Fermi-contact 
one. The Knight shift is dominated by the Fermi contact 
polarization both on B and on Mg. After these calcula- 
tions were completed, we learned about similar calcula- 
tions for the valence electrons relaxation rate on B from 
the Ames groupE3, with the results consistent with those 
reported here. 
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